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ABSTRACT 

With the ASCA X-ray satellite, we perform spatial and spectral studies of the shell-type supernova rem- 
nant (SNR) 7 Cygni that is associated with the brightest EGRET unidentified source 3EG J2020+4017. 
At energies below 3 keV the bulk of the X-ray flux from the remnant is well described by the emis- 
sion of thermal plasma with characteristic temperature kT e ~ 0.5-0.9 keV. In addition to this thermal 
emission, we found an extremely hard X-ray component from several clumps localized in the northern 
part of the remnant. This component, which dominates the X-ray emission from the remnant above 
4 keV, is described by a power-law with a photon index of T ~ 0.8-1.5. Both the absolute flux and the 
spectral shape of the nonthermal X-rays cannot be explained by the synchrotron or inverse-Compton 
mechanisms. We argue that the unusually hard X-ray spectrum can be naturally interpreted in terms 
of nonthermal bremsstrahlung from Coulomb-loss-flattened electron distribution in dense environs with 
the gas density about 10 to 100 cm~ 3 . For given spectrum of the electron population, the ratio of the 
bremsstrahlung X- and 7-ray fluxes depends on the position of the "Coulomb break" in the electron 
spectrum. Formally, the entire high energy 7-ray flux detected by EGRET from 7 Cygni could originate 
in the hard X-ray regions. However, it is more likely that the bulk of 7-rays detected by EGRET come 
from the radio-bright and X-ray dim cloud at southeast, where the very dense gas and strong magnetic 
field would illuminate the cloud in the radio and 7-ray bands, but suppress the bremsstrahlung X-ray 
emission due to the shift of the "Coulomb break" in the electron spectrum towards higher energies. 



Subject headings: acceleration of particles — cosmic rays- 
waves — supernova remnants 



radiation mechanisms : non-thermal — shock 



1. INTRODUCTION 

The origin of the unidentified EGRET 7-ray sources in 
the Galactic plane has been puzzling since their discovery. 
The third EGRET catalog (Hartman et al. 1999) lists 57 
unidentified sources at low Galactic latitude (|6| < 10°). 
Rotation-powered pulsars are likely to account for part of 
these GeV 7-ray sources; five pulsars have been detected 
at GeV energies to date. Another possible origin for some 
of the unidentified sources is the emission from accelerated 
cosmic-rays at the shock of shell-type supernova remnants 
(SNRs). It is reported that the probability to find EGRET 
unidentified sources in the vicinity of shell-type SNRs is 
significantly high (Sturner & Dermer 1995). Although 
the relatively young SNRs, being in their Sedov phase, 
are natural sites of high energy 7-ray production through 
electron bremsstrahlung and hadronic interactions, it has 
been recognized that in most cases the expected 7-ray 
fluxes at MeV/GeV energies are too low to be detected 
by EGRET (Drury, Aharonian, & Volk 1994). However, 
the 7-ray fluxes can be dramatically enhanced in SNRs 
having dense gas environments, e.g., in large molecular 
clouds overtaken by supernova shells (Aharonian, Drury, 

6 Volk 1994). Remarkably, among the SNRs possibly de- 
tected by EGRET are the radio-bright and nearby objects 

7 Cygni, IC 443, W44, and W28 (Esposito et al. 1996) 
that are all associated with molecular clouds. The recent 
searches for SNRs in the vicinity of some of the unidenti- 
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fled EGRET sources revealed new evidence of 7-ray emis- 
sion from "supernova remnant-molecular cloud" interact- 
ing systems (Combi, Romero, & Benaglia 1998; Combi et 
al. 2001). If gamma-rays of such systems have hadronic 
origin, and the acceleration of protons extends to 10 TeV 
and beyond, the TeV 7-ray emission from these objects 
should be detectable also by current ground-based detec- 
tors (Aharonian, Drury, & Volk 1994). At energies below 
the threshold of production of 7r°-decay 7-rays around 70 
MeV, the electron bremsstrahlung remains the only notice- 
able 7-ray production process. Therefore, the detection of 
7-rays by EGRET down to several tens of MeV implies the 
existence of low-energy (< 100 MeV) electrons in these ob- 
jects. The best information about higher energy electrons, 
typically between 1 and 10 GeV electrons, is provided by 
the synchrotron radio emission at GHz frequencies. 

The accelerated electrons produce also nonthermal 
X-rays through the synchrotron radiation and the 
bremsstrahlung. Both channels contain unique informa- 
tion about the nonthermal electron populations at ex- 
tremely high (multi-TeV), and very low (sub- MeV) ener- 
gies, respectively. Thus, the X-ray observations may help 
to reveal the nature of the unidentified EGRET sources 
that are associated with SNRs by looking for ultrarela- 
tivistic and subrelativistic electrons, and probing the en- 
vironments of the remnants. 

The 7 Cygni (G78. 2+2.1) SNR has a clear position- 
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correlation with the brightest unidentified 7-ray source 
3EG J2020+4017. It is a nearby (1-2 kpc) shell-type 
SNR with the radio shell of ~ 60' diameter (Higgs et 
al. 1977). The radio flux density of 340 Jy at 1 GHz 
ranks it as the fourth brightest SNR in the sky at this 
frequency (Green 2000). Almost 60% of the radio flux 
comes from southeastern part which has been known as 
DR4 (Downes & Rinehart 1966). The spectral index 
of the radio spectrum averaged over the whole remnant 
is measured as a ~ 0.5 (Green 2000). Its variation 
across the remnant is as small as Aa ~ ±0.15 (Zhang 
et al. 1997). The EGRET source has the steady flux of 
F(E > 100 MeV) = (12.6 ±0.7) x 10~ 7 photons cm^s" 1 
and a best-fit power-law index of 2.07 ± 0.05 (Esposito 
et al. 1996). Prior to the EGRET detection, 7-ray source 
2CG078+2 was detected in the vicinity of 7 Cygni with the 
COS-B satellite (Pollock 1985). Despite extensive searches 
for TeV 7-ray emission, significant excesses have not been 
detected (Buckley et al. 1998). Since a simple extrapola- 
tion of the EGRET flux exceeds the Whipple upper limit 
by an order-of-magnitude, the spectrum must have a cutoff 
or steepen well below the TeV energy (Gaisser, Protheroe, 
& Stanev 1998; Buckley et al. 1998). 

Presented here are the results from and implications of 
ASCA observations of the 7 Cygni SNR. Of particular in- 
terest is the discovery of hard X-ray emission localized to 
the northern part of 7 Cygni. The plan of this paper is as 
follows. The ASCA observations are briefly summarized 
in §2. We perform the X-ray image analysis in §3.1 and 
the spectral fits in §3.2. In §4.2, we interpret the soft X- 
ray data in terms of shock-cloud interaction. The origin 
of the hard X-ray emission is discussed in §4.3. 

2. OBSERVATIONS 

We performed ASCA observations of 7 Cygni SNR twice 
in 1996 and 1997. The northern part of the remnant was 
observed for 40 ks in 1996 May. Observations were carried 
out in 1997 May with three pointings toward the north, 
south, and east of 7 Cygni with 60, 16, and 12 ks duration, 
respectively. The data from two Gas Imaging Spectrome- 
ter (GIS; Makishima et al. 1996) detectors were acquired 
in the standard pulsc-hcight mode. Two Solid-state Imag- 
ing Spectrometer (SIS; Burke et al. 1994) detectors were 
operated in 4-CCD mode. Each detector is coupled to its 
own telescope with nested conical foil mirrors. Once com- 
bined all observations, field of view (FOV) of GIS covers 
almost all part of the remnant. 

The data from the detectors were screened by standard 
procedures, except for the strict rise-time screening in the 
GIS analysis, because of the mode we chose for the opera- 
tion. The observed count rate for the north of 7 Cygni was 
0.57 (0.59) counts/s within a central 20' radius in the 0.7- 
10 keV energy band with GIS2 (GIS3). The degradation 
of the SIS performance due to the accumulated radiation 
damage in orbit is found to be substantial in our data. 
Even for the data taken with faint mode, in which RDD 
(Residual Dark Distribution) correction can be applied, 
the spectral resolution of SIS is found to be of no advan- 
tage to that of GIS. In addition, the smaller effective area 
at high energies and the narrower FOV for SIS than those 
for GIS do not serve our purpose. The results presented 
in this paper are from the GIS data. 



3. DATA ANALYSIS AND RESULTS 

3.1. X-ray Images 

The ASCA GIS data for the four pointings were com- 
bined to construct the X-ray images in selected energy 
bands. We discard the data outside the central region of 
the detector with a radius greater than 20'. By utilizing 
the night-earth observation data set, we subtract the in- 
strumental background from each pointing data. Since we 
did not obtain the rise-time information, the normalization 
of the background image is increased by 30% as compared 
with the nominal case. The background-subtracted images 
are combined with corrections for the exposure time and 
the vignetting effects. Finally the images are smoothed 
with a function comprised of a narrow core (< 1') and a 
broad wing to simulate the point spread function (PSF) 
of the telescope. Figure la, 16, and lc show the resultant 
X-ray images in 0.7-1 keV (soft), 1-3 keV (medium), and 
4-7 keV (hard) energy bands, respectively. 

We find strikingly different appearance of the SNR in 
these energy bands. In the soft band image, clumpy 
sources appear at north and limb-brightened emission is 
seen at south. The medium image shows widespread emis- 
sion distributed from north to southeast, in addition to the 
southern bright shell. In the hard energy band, several 
clumpy sources in the north region stand out dramatically 
(Figure lc). Hereafter, we refer to the clump centered at 
a celestial coordinate (J2000) of R.A.= 20 h 21 m 12 s , decl.= 
40°47'55" as CI and the other one at R.A.= 20 h 20 m 00 s , 
decl.= 40°40'03" as C2. Since another hard source C3 is 
located at the western edge of the FOV, we do not report 
a spectral study of C3 in this paper. 

The hard source C2 appears to be extended by a few 
arcminutes. In order to evaluate quantitatively, we com- 
pared the count rate profile centered on the clump C2 with 
the detector PSF in the 3-10 keV band, as plotted in Fig- 
ure 2. The background level shown there is taken from the 
southern part of the remnant. It is obvious that the hard 
X-ray emission extends up to 4' -6' beyond the radial pro- 
file of a single point source. In Figure 2, we also plot the 
simulated profile of a circular source which has uniform 
surface brightness with a radius of 4' centered on C2. On 
the other hand, we obtain no sign of spatial extension of 
the source CI. 

In order to compare the X-ray distribution with radio 
synchrotron emission, the X-ray images on gray scale are 
overlaid with the 4850 MHz radio contours in Figure 1. 
The X-ray sources at the north coincide with the radio- 
bright region and the southern X-ray shell with a fainter 
radio arc. The X-ray morphology in general bears close 
resemblance to the radio on large scales, but the X-ray 
emission in the proximity of the region R2 appears to an- 
ticorrelate with the radio emission. The DR4, brightest 
radio region at the southeast, is quite dim in the X-ray 
band. It is noticeable that there are enhancements of ra- 
dio emission possibly associated with the hard sources CI 
and C2. 

The second (Thompson et al. 1995) and third (Hart- 
man et al. 1999) EGRET catalog report somewhat dif- 
ferent positions for the unidentified 7-ray source inside 
7 Cygni. In Figure lc we superimpose their 95% confi- 
dence error contours which are approximated as circles. 
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Whereas Brazier et al. (1996) reported a point-like X- 
ray source RX J2020. 2+4026 by ROSAT close to the 
remnant center and within the EGRET error circle, no 
point sources were detected there by our ASCA obser- 
vations. The flux of RX J2020.2+4026, estimated to be 
~4x 10~ 14 ergs cm _2 s _1 assuming a photon index of 2, 
is indeed below the sensitivity of our observation. In §4.3 
we discuss a possible connection between the hard X-ray 
sources and the EGRET unidentified 7-ray emission. 

3.2. X-ray Spectra 

For the purpose of spectroscopic studies of the X-ray 
emissions associated with the 7 Cygni SNR, contaminat- 
ing emission that is overlapping the remnant must be taken 
into account adequately because the line of sight toward 
7 Cygni passes along the Orion-Cygnus spiral arm. Be- 
cause of the lack of a blank field in our observations, we 
analyzed four ASCA archive data in neighboring fields 
(Field 1-4 in Table 1) to estimate the X-ray emission at 
the 7 Cygni field unrelated to the SNR itself. For each 
field, the GIS spectrum is obtained by integrated over the 
central detector region with a radius of 20' after eliminat- 
ing resolved sources. Each spectrum is subtracted by the 
high-latitude blank-sky spectrum as a sum of the cosmic 
X-ray background and the instrumental background. Fig- 
ure 3 shows the resulting spectra in the 1.2-2.5 keV and 
3.5-8 keV band, after the correction for the integration sky 
area. We find similar energy spectra from these neighbor- 
hood fields, except for Field 2 where the Galactic latitude 
is highest and the Galactic column density is lowest as 
compared with the other fields. Once Field 2 is excluded, 
field-to-field variations of spectral data in the soft band are 
found to be insignificant. The hard 3.5-8 keV spectrum of 
the southern field of 7 Cygni is in agreement with those of 
Field 3 and 4 within their statistical uncertainties. Conse- 
quently, we consider that Field 3 and 4 provide us a good 
approximation of the contaminating emission that should 
be subtracted from the spectral data of 7 Cygni. Since no 
bright sources are found in the ASCA data of Field 4, we 
have chosen Field 4 as "background" field. 

We derive energy spectra of regions R1-R3 and clumps 
CI and C2. The clump spectra are extracted from the cir- 
cular regions of a radius of 6'. We exclude photons falling 
within 6' radius centered on CI and C2 from the spectrum 
of the region Rl. Each accumulated on-source spectrum 
is subtracted by the "background" Field 4 spectrum that 
is extracted from an identical detector region to each on- 
source data. To improve statistics, spectra of two GIS de- 
tectors are always added. Background-subtracted spectra 
of R1/R3 and C2 are shown in Figure 4 and Figure 5, re- 
spectively. Several emission lines of Mg K (~ 1.4 keV) and 
Si K (~ 1.9 keV) are evident in every spectra, indicative 
of thin thermal plasma with a typical temperature of « 1 
keV. Remarkably, the spectrum of C2 exhibits very flat 
continuum emission above 3 keV. There are known errors 
in the calibration of the conversion of the photon energies 
to the pulse invariant channels of the GIS detector below 
the xenon-L edge of 4.8 keV. In the following spectral fit- 
ting, an artificial energy shift of —50 eV to each applied 
model is introduced to alleviate these errors (Buote 1999). 

We have first attempted to fit the 0.7-8 keV spectrum 
of R3 by a thin thermal plasma model (Mewe, Gronen- 



schild, & Oord 1985; Liedahl et al. 1990) in which col- 
lisional equilibrium ionization (CEI) is realized. Photo- 
electric absorption along the line of sight is taken into ac- 
count using the cross-sections from Morrison & McCam- 
mon (1983). Elemental abundances are fixed to the solar 
values of Anders & Grevesse (1989) throughout this pa- 
per unless otherwise mentioned. The CEI plasma model 
cannot give an acceptable fit owing to a large discrepancy 
between the actual data and the model between Mg and 
Si K emission lines. Even when the abundances of alpha 
elements Ne, Mg, Si, S, and Fe are individually allowed to 
vary ranging from 0.1 to 10 solar, the spectral data can- 
not be described properly. In order to model emission line 
features, we take account of the effects of non-equilibrium 
ionization (NEI; Itoh 1979) by adopting a plasma emission 
code based on Masai (1994). In an NEI plasma, degree 
of ionization and consequently line emissivities depend on 
the ionization timescale n e t, where n e represents the elec- 
tron density and t the passage time after shocked. The 
NEI plasma model yields an acceptable fit with a reduced 
X 2 (d.o.f) = 1.19(26). The best-fit temperature and ion- 
ization timescale with la errors are kT e — 0.76^Jg keV 

and n e t = 5.81+4 x 10 10 cm~ 3 s, respectively. 

We found the Rl spectrum, as compared with R3, shows 
distinctive features, namely a strong emission at ~ 0.9 keV 
and a hard continuum above 3 keV. The former is consis- 
tent with the fact that Rl is very bright particularly in 
the 0.7-1 keV energy band; the latter could be contam- 
ination by the hard sources CI and C2. Regarding the 
spectral fit of Rl , we include helium- like neon (Ne IX) Ka 
line (0.923 keV), in addition to the CEI plasma model that 
predominantly describes the 1-3 keV emission. The CEI 
plasma plus Ne IX line model, however, cannot give a sta- 
tistically acceptable fit, owing to residuals by the hard con- 
tinuum. Then, by adding power-law as a third component 
to describe the hard continuum, we obtain an acceptable 
fit for the Rl spectrum as summarized in Table 2. The 
best-fit temperature is kT e — 0.56^^05 keV; tne photon 
index of the power-law component is T = 1.2 Ijjj'g. 

The spectrum of the region R2 is fitted by a CEI plasma 
model alone. The simple model yields an acceptable fit 
with the temperature of kT e = 0.53 ± 0.07 keV, which 
is in complete agreement with the value obtained for the 
region Rl. 

The spectral data of the clump C2 indicate the presence 
of a hard continuum emission in the 3-8 keV energy band. 
The hard emission is seen in the spectrum of CI as well. As 
the case of Rl, we employ a three-component model com- 
prised of the Ne IX line (0.92 keV), CEI thermal plasma, 
and power-law spectrum. We find that a good fit cannot 
be obtained if we omit the power-law component. We fit- 
ted the three-component model to the 0.7-8 keV spectral 
data of the clumps CI and C2, by freezing the plasma tem- 
perature to the best-fit value of the region Rl and atten- 
uating all components with a common absorption column 
of A H = 0.84 x 10 22 cm~ 2 that is also obtained for Rl. 
The power-law components are found to be very flat; the 
best-fit photon indices are T = 1.5 ± 0.5 and 0.8 ± 0.4 for 
CI and C2, respectively. Instead of adding a power-law, 
we also attempted to add a thermal bremsstrahlung com- 
ponent to model the high-energy part of the C2 spectrum. 
The 90% lower-limit on the temperature of the thermal 
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bremsstrahlung emission is set to be kT e = 5.8 keV. 

4. DISCUSSION 

4.1. Estimation of SNR Distance and Age 

Higgs et al. (1977) derived the distance to 7 Cygni as 
1.8 ± 0.5 kpc based on the S-D relation which is a sta- 
tistical property of the radio brightness of SNRs. Lan- 
decker, Roger & Higgs (1980) estimated as 1.5 ± 0.5 kpc 
and pointed out the progenitor of the 7 Cygni remnant 
was possibly a member of the Cyg OB9 association at 
1.2 ±0.3 kpc. The absorbing column density provides addi- 
tional information about the distance. Maeda et al. (1999) 
reported Wolf-Rayet binary V444 Cyg located close to the 
7 Cygni sky field is attenuated by the interstellar column 
density of 7V H = (1.1 ± 0.2) x 10 22 cm" 2 , similar to the 
7 Cygni remnant. Thus, we expect the distance to 7 Cygni 
is probably close to the distance to V444 Cyg, 1.7 kpc. In 
view of these arguments, we take the distance D = 1.5 kpc 
as most probable value. 

The X-ray distribution of the region R3 shows clear 
arc-like morphology along the outer boundary of the shell 
structure observed at radio frequencies. The X-ray spec- 
trum was modeled by a single temperature thermal emis- 
sion with the solar abundance. These features are indica- 
tive of a thin thermal plasma in the immediate postshock 
region of the primary blast wave propagating through the 
interstellar matter. 

Assuming the equipartition between the shocked elec- 
trons and ions and Rankine-Hugoniot jump conditions for 
a strong shock, the electron temperature is related to the 
shock velocity v s as kT e = (3/16)/!to p i> 2 , where /i — 0.6 is 
a mean mass per particle in units of the proton mass m p . 
The best-fit value of kT e = 0.761q'qq keV corresponds to 
the shock velocity v s = 800lg° km/s. 

Provided that the remnant is in the Sedov adiabatic ex- 
pansion phase, the age of 7 Cygni can be estimated to be 
T ago = (2/5)R/v s , where R is the radius of the supernova 
shock front. Given the angular radius of 9 ~ 30', we have 
the physical radius R = DO c± 13.5Z?i.5 pc, where D1.5 is 
the distance in units of 1.5 kpc. Thus, we obtain an age 
estimate of T ago — 66OOD1.5 yr. 

4.2. Signatures of Shock-Cloud Interaction 

On the basis of the H I line emission and absorption to- 
ward the 7 Cygni SNR, Landecker, Roger & Higgs (1980) 
suggested that the supernova explosion took place in the 
slab of an interstellar cloud oriented north to southeast. 
The soft (1-3 keV) X-ray emission "belt" from north to 
southeast appears to agree fairly well with the spatial dis- 
tribution of the H I line features. A possible scenario ex- 
plaining this coincidence is that the X-ray emission belt is 
due to thermal evaporation of clouds (White & Long 1991) 
as a consequence of the collision between the supernova 
blast wave and the dense clouds. The supernova shock 
may expand into a cavity produced by progenitor stellar 
winds, and then encounter the cavity wall comprised of 
the clouds (Chevalier 1999). 

We found the possible evidence of the strong emission 
line of helium-like Nc from the Rl spectrum. The X-ray 
map in the 0.7-1 keV band suggests clumpy nature of the 
Ne IX line emission. The upper limit for the hydrogen-like 



Ne line intensity is found to be 6.7% of the helium- like 
intensity. This gives an upper limit of 0.24 keV for the 
temperature of plasma emitting the Ne IX line if we as- 
sume ionization equilibrium state. The low-temperature 
and clumpy appearance are indicative of the shock-cloud 
interactions at the northern part of the remnant. 

The hard sources CI and C2 reside in the northern re- 
gion where the shock-cloud interaction is likely to occur. 
Therefore, we identify the hard sources with shocked dense 
cloudlets. This speculation may be encouraged by the 
presence of high-velocity H I cloudlets with a density of 
10-100 cm~ 3 and an apparent radius of 4'-6' in the south- 
east of the remnant (Landecker, Roger & Higgs 1980). Let 
us calculate the density of plasma that is responsible for 
the derived Ne flux /nc ~ 6(0.8) x 10~ 3 photons cm _2 s _1 
of the clump C1(C2) with reasonable parameters. Given a 
temperature of 0.1 keV and a spherical volume of an angu- 
lar radius of 2' (assuming distance D = 1.5 kpc), we have 
the density n 60(8) cm~ 3 by using the Ne line emissivity 
of Mewe, Gronenschild, & Oord (1985), in agreement with 
the typical density of H I cloudlets. 

4.3. On the Origin of Hard X-ray Emission 

An important finding of our observations is the clumpy 
hard X-ray emissions from the regions CI and C2 (here- 
after, HXC - the hard X-ray clump). We found that C2 
has an extended structure with an effective radius of about 
4'. The HXC components are extremely hard and charac- 
terized by a power-law with T ~ 0.8-1.5. As we argued 
above, it is likely the HXC originates in shocked dense 
cloudlets with a gas density between 10 and 100 cm~ 3 . 

The thermal-bremsstrahlung interpretation of the high- 
energy part of the C2 spectrum requires very high temper- 
ature (kT e > 5.8 keV). Since the shocks of an evolved SNR 
are not sufficiently energetic to heat the gas to such tem- 
peratures, the nonthermal origin of the hard X-ray emis- 
sion seems a more favorable option. The synchrotron radi- 
ation, the inverse-Compton scattering, and the nonthermal 
bremsstrahlung of shock-accelerated electrons are three 
natural production mechanisms of hard X-rays. The over- 
all flux of the HXC in the 2-10 keV interval is estimated 
as Fx — 4.5 x 10 _12 ergs cm _2 s _1 , which corresponds to 
the source luminosity Lx — 1-2 x 10 33 Z? 2 5 ergs s _1 . 

The nonthermal luminosity of this SNR peaks at high 
energy 7-rays, if the association with the EGRET source is 
real. The 7-ray spectrum measured by EGRET is shown 
in Figure 6, together with the HXC spectrum. The ob- 
served photon flux of the EGRET source is translated into 
a luminosity of L 7 ~ 1.4 x 10 35 L> 2 5 ergs s" 1 (100 MeV- 
2 GeV) for the photon index 2.1. The reported posi- 
tions of the EGRET 7-ray source 2EG J2020+4026 and/or 
3EG J2020+4017 do not coincide exactly with the HXC. 
However the large systematic errors in the EGRET posi- 
tion, which strongly depend on the chosen diffuse 7-ray 
emission model (Hunter et al. 1997), do not allow certain 
conclusions concerning the location of the 7-ray produc- 
tion region. Therefore it is interesting to test whether the 
HXC can be a counterpart of the unidentified GeV 7-rays 
source on theoretical ground. 

4.3.1. Coulomb Cooling and Nonthermal Bremsstrahlung 

The emission region of the HXC is supposed to be the 
dense cloudlet engulfed by the supernova blast wave, and 
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therefore the density is probably high (n <~ 10-100 cm~ 3 ). 
Below we assume that accelerated electrons are contin- 
uously injected in the HXC, by the accelerator(s) being 
located cither outside or inside the HXC. Possible internal 
accelerator could be the secondary shocks generated by the 
interaction between the primary blast wave and the dense 
cloudlets. Indeed, Bykov et al. (2000) considered elec- 
tron acceleration at a slow shock (the order of 100 km/s) 
transmitted inside a dense cloudlet. A tail shock could be 
formed behind the cloudlets (Jones & Kang 1993) as well, 
although an efficient particle acceleration at the tail shock 
is not obvious. On the other hand, the external accelerator 
should be the primary shock of the SNR. 

The energy spectrum of accelerated electrons depends 
on the acceleration mechanism, the discussion of which is 
beyond the scope of this paper. For simplicity we assume 
the electron production rate in momentum space is de- 
scribed by a single power-law function Q(p) oc jT K . Such 
a function implies power-law distributions of electrons also 
in the kinetic energy space, although in the non-relativistic 
and ultrarelativistic regimes the power-law indices are dif- 
ferent, namely Q(E) oc £-( s +!)/ 2 a t E < m e c 2 while 
Q(E) oc E~ s at E ^> m e c 2 , where s = k + 2. In the case 
of diffusive shock acceleration at the shock front, accord- 
ing to the standard model (Bell 1978; Blandford & Eichler 
1987), the spectral index is determined by the compression 
ratio of the shock as s = (r + 2)/(r — 1). For strong shocks 
with r ~ 4 the spectral index is close to 2. 

For simplicity we suppose that the electrons are effec- 
tively trapped in HXC; thus the present-day electron spec- 
trum is formed by accumulation of freshly accelerated (or 
arriving from external accelerators) electrons during the 
entire history of the HXC. If the energy-losses can be ne- 
glected, the electron spectrum becomes N{E) — r agc Q(E), 
where the age of the accelerator is approximated by the 



age of 7 Cygni, i.e. 



7000 yr. It is valid, how- 



ever, only in the intermediate energy region, typically be- 
tween 300 MeV to 10 GeV, where the bremsstrahlung 
dominates the radiative losses. Indeed, the characteris- 
tic lifetime of electrons against the bremsstrahlung losses, 
Tbr = E/(-dE/dt) hr ~ 4.3 x 10 7 (n/10 cm" 3 )- 1 yr is con- 
siderably larger than the age of 7 Cygni as long as the gas 
density does not exceed 6 x 10 4 cm -3 . 

The high density and presumably strong magnetic field 
in the HXC make, however, the cooling processes cru- 
cial in the formation of the electron spectrum at low and 
very high energies. The relevant cooling processes are the 
Coulomb (or, in the neutral gas, ionization) losses for low- 
energy electrons and the synchrotron losses for high-energy 
electrons. The inverse-Compton losses could be omitted if 
the magnetic field in HXC exceeds 10 /LtG. Also, it should 
be noticed that even in case of significant bremsstrahlung 
losses (Tbr < T a g o), the process does not modify the elec- 
tron spectrum. 

Below several hundreds MeV, the energy losses of elec- 
trons are dominated (independent of the gas density) by 
the Coulomb losses (see e.g. Hayakawa 1969). The cool- 
ing time of electrons due to the Coulomb interactions is 
given by (Rephaeli 1979) 



E 



^cou — 



(3/2)n e ccr T (TOeC 2 )/3- 1 lnA 



(1) 



4.3 x 10' 



iJv) yr ' 



(2) 



,10 cm- 3 , 

where ctt is the Thomson cross-section, c is the velocity 
of light, E and (3 are the kinetic energy and the velocity 
of nonthermal electrons in units of c, m e is the electron 
mass, and In A is the Coulomb logarithm which is set to 
be 40. Equating t cou and T age gives the break energy E cou , 
"Coulomb break", below which the electron spectrum is 
flattened: 



E n 



1.6 



MeV, 



(3) 



.10 cm- 3 J V7000 yr, 
where we use the approximation of /? ~ 1. At energies 
below E cou the Coulomb losses significantly modify the 
acceleration spectrum, N(E) = t cou (E)Q(E); namely, at 
relativistic energies, N(E) oc E~ s+1 . 

The differential energy spectrum of the bremsstrahlung 
emission from these electrons is calculated as (see e.g. Blu- 
menthal & Gould 1970) 

1(e) ^ J dEN(E)cf3 (n H ^f + n Hc ^f + n e ^ ' 

6 (4) 

where nn, nu e , and n e are hydrogen, helium, and elec- 
tron number densities, respectively, and da/de is the dif- 
ferential cross-section for emitting a bremsstrahlung pho- 
ton in the energy interval e to e + de. We adopt ratios 
nne/n-ft = 0.1 and n e /nn = 1.2. In the ultrarelativis- 
tic regime the electron-electron bremsstrahlung becomes 
comparable to the electron-proton bremsstrahlung. Since 
the bremsstrahlung cross-section is in inverse proportion to 
the emitted photon energy, da/de oc and only slightly 
(logarithmically) depends on the electron energy in the 
relativistic regime, the bremsstrahlung photon spectrum 
almost repeats the power-law spectrum of parent elec- 
tron distribution of N(E) = kE~ a , i.e. 1(e) oc e~ a if 
a > 1. Then, we should expect a broken power-law spec- 
trum with 1(e) oc e~ s at e > E cou , and 1(e) oc e~ s+1 at 
e < E cou in the relativistic domain. The energy spectrum 
of bremsstrahlung photons is essentially determined by the 
gas density and the age. Below the "Coulomb break" , for 
any reasonable acceleration index of s < 2.5, we see a 
hard differential spectrum with a photon index less than 
1.5. It should be stressed that, for the electron spectrum 
harder than E^ 1 (E~ x l 2 in the non-relativistic case), the 
bremsstrahlung photons obey a standard e^ 1 type spec- 
trum. In the ASCA energy band, we should always expect 
the standard spectrum. 

4.3.2. Synchrotron Radiation 

The synchrotron cooling time of electrons in the mag- 
netic field B is given by 



7TO e C 



syn (4/3)ca T C/ B 7 2 
B 



= 1.2 x 10 3 



10- 5 G 



10 14 eV 



(5) 



yr, (6) 



where 7 is the electron Lorentz factor and Ub = B 2 / (8ir) 
is the energy density of the magnetic field. Equating r syn 
and T agc gives the break energy E syn above which the elec- 
tron spectrum is steepened: 



18 



B 



'ago 



io- 5 g; V 700 °y r 



TeV. (7) 
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Thus if the maximum energy of electrons exceeds E syn , 
we should expect a break in the electron spectrum. In 
particular, in the case of power-law acceleration spec- 
trum, the loss-steepened spectrum has a power-law form 
N(E) =r syn {E)Q{E) (xB— 1 . 

The average synchrotron photon energy radiated by an 
electron with energy of En (in 10 14 eV) in the field B- 5 
(in 1(T 5 G) is e ~ 1.6-B_ 5 -E 2 4 keV. Then, the break en- 
ergy E syn at which the synchrotron lifetime is balanced 
with the SNR age is translated into the synchrotron pho- 
ton energy of e syn - 50 (-B/10- 5 G)- 3 (T age /7000 yr)~ 2 eV, 
above which noticeable spectral steepening should appear. 
Therefore the observed flat X-ray spectrum up to 10 keV 
would require unrealistically small magnetic field in the 
HXC region, B < 2 ^G. 

Furthermore, it is easy to show that in the case of ac- 
celeration of electrons inside HXC, independent of the 
strength of the magnetic field, the cut-off energy in the 
synchrotron spectrum cannot exceed a few eV. Indeed, if 
the maximum energy of accelerating electrons is limited 
by the synchrotron losses, the balance between the accel- 
eration and cooling timescales determines the maximum 
electron energy. The timescale of the first-order Fermi ac- 
celeration in the Bohm diffusion limit is given by 

Tacc ~ — o- (8) 

where tl = m e j/(eB) is the Larmor radius of elec- 
trons and v s is the shock velocity. Equating r syn in 
equation (5) and T acc gives the maximum energy E m ~ 
6.2 ( J B/10~ 5 G)~ 1 / 2 (w s /100 km/s) TeV, and its character- 
istic photon energy is e m ~ 3.8 (v s /100 km/s) 2 eV that is 
independent of the magnetic field strength. Hence, regard- 
less of the magnetic field strength, the first-order Fermi 
acceleration mechanism cannot explain the extremely flat 
X-ray spectrum by synchrotron radiation, unless the shock 
speed exceeds 5000 km/s. The latter condition obviously 
cannot be satisfied in 7 Cygni. 

4.3.3. Inverse- Compton Emission 

The hard X-rays could be produced also by inverse- 
Compton (IC) scattering of seed photons off relativistic 
electrons. In the 7 Cygni SNR, the cosmic microwave 
background (CMB) radiation dominates over other ambi- 
ent photon fields like the diffuse Galactic infrared/optical 
radiation and the infrared emission from shock-heated dust 
(Gaisser, Protheroe, & Stanev 1998). 

Relativistic electrons which boost the CMB photons up 
to X-rays also emit synchrotron radiation in the radio 
band. Characteristic energies of synchrotron photon e syn 
and IC photon ei c (in keV) produced by the same electrons 
are related as (e.g. Aharonian, Atoyan, & Kifune 1997) 

e syn ~ 0.7 x 10~ 7 e icJ B_ 5 eV (10) 

for the CMB seed field. The ratio of the spectral power 
/(e) = e 2 x 1(e) of the IC X-ray emission to the syn- 
chrotron emission at the relevant radio frequency is 

/ic(eic) _ U CMB 
/syn(e S yn) Ub 



where C/cmb = 0.25 eV cm -3 is the energy density of the 
CMB radiation. The spectral power of the observed radio 
emission from the whole remnant, in the frequency range 
from 408 MHz to 4.8 GHz (Wendker, Higgs, & Landecker 
1991; Higgs et al. 1977) is well described as 

/radio(e) - 5 X 10- 12 ( 1Q -5 c y ) " ^ CT S S ™ 2 S^. (12) 

Since almost 60% of the radio flux comes from the 
southeast of 7 Cygni, from equation (10) and (11) we 
find that the IC spectral power cannot exceed 1.7 x 
10 _14 e-' c 5 -BZ5' 5 ergs cm _2 s _1 . This upper-limit is less 
than 1% of the flux of the hard X-ray emission detected. 
Thus, the IC radiation hardly can explain the flux of HXC, 
unless the magnetic field is extremely weak (< 1 fiG). 

4.3.4. Nonthermal Bremsstrahlung V/7 -radiation 

The nonthermal bremsstrahlung from the accelerated 
electrons is a natural source of the HXC flux, because the 
shocked dense cloudlets act as an effective target for ener- 
getic electrons. Because of Coulomb interactions the high 
density gas leads to significant hardening of low-energy 
electrons below the "Coulomb break" , giving rise to the 
standard e _1 type bremsstrahlung spectrum at the X-ray 
band, which perfectly agrees with the ASCA data. The 
bremsstrahlung spectrum above the "Coulomb break" es- 
sentially repeats the acceleration spectrum of electrons. 
For the given spectral index of accelerated electrons, the 
ratio of the X- and 7-ray fluxes depends only on the 
"Coulomb break" energy which in its turn depends on the 
gas density and the age of the accelerator. 

In Figure 6 we show the results of numerical calcula- 
tions for two sets of parameters which describe the gas 
density and the acceleration spectrum of electrons, by as- 
suming the electron bremsstrahlung is responsible for both 
the ASCA hard X-ray and the EGRET 7-ray fluxes. For 
the electron spectrum with the acceleration index s = 2.1, 
the best fit is achieved for a gas density of n = 34 cm -3 . A 
steeper acceleration spectrum with s — 2.3 requires larger 
gas density, n — 130 cm~ 3 . Note that the adopted ac- 
celeration spectra are consistent with the reported radio 
spectral index a — 0.5 ± 0.15. We also suppose an expo- 
nential cutoff in the electron spectrum at 10 TeV. 

If the electron distribution with s = 2.3(2.1) extends be- 
yond GeV energies, for the magnetic field 10 -5 G the cal- 
culated radio flux density amounts to about 10% (60%) of 
the measured radio flux density integrated over the whole 
remnant. Furthermore, if the electron distribution extends 
beyond TeV, we found the bremsstrahlung spectrum with 
,s = 2.1 exceeds the Whipple upper-limit, whereas the 
spectral index of s = 2.3 is still in agreement with the 
Whipple data. Meanwhile both combinations of model pa- 
rameters satisfactorily fit the spectral shape and the abso- 
lute flux of hard X-rays. The ignorance of energy losses of 
electrons would lead to significantly steeper X-ray spectra, 
and would also result in overproduction of absolute X-ray 
fluxes (dotted curves in Figure 6). Note that the main 
contribution to X-rays comes from relatively high energy 
electrons with energies close to 10 MeV for the electron 
index s — 2.3. 

Because of poor angular resolution, the EGRET mea- 
surements do not provide a clear information about the 
site(s) of production of high energy 7-rays. Nevertheless, 
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it is likely that only a part (perhaps, even only a small 
part) of the reported high energy 7-ray fluxes originates in 
the HXC region. The 7-ray fluxes could be suppressed by 
assuming lower gas densities. Indeed, such an assumption 
would lead to the shift of the "Coulomb break" energy in 
the electron spectrum to lower energies, and the predicted 
high energy 7-ray spectra would appear significantly below 
the reported EGRET fluxes (a solid curve in Figure 6). 

A more likely candidate for production of the bulk of 
high energy 7-rays is the region called DR4 from which 
most of the radio emission emerges. A massive cloud with 
a density of ~ 300 cm~ 3 occupying <~ 5% of the SNR vol- 
ume has been suggested to exist in the vicinity of DR4 
to explain the 7-ray flux (Pollock 1985). Actually the 
EGRET error circle reported is somewhat away from the 
HXC but closer to the DR4. A density of - 300 cm" 3 
is higher than the upper limit density of the HXC. Such 
high gas density implies a high (about 50 McV) "Coulomb 
break" energy in the electron spectrum, and therefore con- 
siderable suppression of X-ray flux. This could naturally 
explain the lack of noticeable hard X-ray fluxes from the 
DR4 region which is bright in radio and possibly 7-rays. 

Finally, we briefly discuss the power consumption due 
to the rapid Coulomb losses. For the gas density of about 
100 cm~ 3 in the HXC, the X-ray flux is produced pre- 
dominantly by electrons with energies of about 10 MeV. 
The X-ray flux is roughly proportional to the product of 
the gas density and the number of relativistic electrons, 
because the relativistic bremsstrahlung cross-section de- 
pends only logarithmically on the electron energy. On 
the other hand, the Coulomb energy loss rate of rela- 
tivistic electrons, dE/dt oc E/t cou , is proportional to the 
gas density and almost independent of the electron en- 
ergy. Therefore the energy loss rate of the bulk of the 
X-ray emitting electrons can be uniquely determined by 
the X-ray luminosity. The measured X-ray luminosity, 
Lx — 1.2 x 10 33 Di 5 ergs s -1 , can be converted to the 
energy loss rate of L e ~ 5 x l0 37 Df 5 ergs s _1 . The energy 
released in relativistic electrons is roughly estimated as 
W e <~ r agc L e ~ 10 49 £>f 5 ergs. This enormous energy de- 
position due to Coulomb collisions would heat the emission 
region of the HXC. Subsequently the heat would be radi- 
ated away in the far-infrared band by molecular line emis- 
sion, if the gas is comprised of molecules. The observed in- 
frared luminosity of 7 Cygni (Saken, Fesen, & Shull 1992) 
is comparable to the energy loss rate estimated above. On 
the other hand, if the emission region of the HXC are pre- 
dominantly shock-ionized plasmas, the deposited energies 
by the accelerated electrons would heat up the plasmas. 

5. SUMMARY AND CONCLUSION 

The X-ray emissions emerging from the 7 Cygni SNR 
are found to be complex. Spatial and spectral studies of 
the X-ray data reveal the presence of four components: 
(1) limb- brightened thermal emission with a temperature 
of ~ 0.8 keV; (2) widespread emission with a temperature 
of ~ 0.6 keV aligned from north to southeast bounded 



by the radio-bright regions; (3) strong emission lines from 
Ne IX ions in the vicinity of the northern bright-radio re- 
gion; (4) clumpy hard emissions, which is best described 
by unusually hard power-law spectral distribution with a 
photon index T ~ 0.8-1.5. 

The limb-brightened component is considered to be 
thermal plasma emission from the immediate postshock re- 
gion of the SNR blast wave that is propagating through the 
interstellar medium. The temperature gives the age esti- 
mate of T ag0 ~ 6600 yr based on the Sedov evolution of the 
remnant. The soft X-ray emissions from north to south- 
east are likely to be caused by the interaction between the 
supernova shock and the cavity wall comprised of ambient 
clouds. In particular, the intense neon line emission would 
be attributable to the low-temperature plasma generated 
by the shock-cloud interaction. 

The extremely hard X-ray emission of the HXC is nat- 
urally explained by the nonthermal bremsstrahlung from 
the loss-flattened electron distribution. By assuming that 
the HXC is a counterpart of the EGRET unidentified 7-ray 
source 3EG J2020+4017, we estimate the density of the 
emission region to be <~ 130 cm~ 3 , for the electron index 
2.3. The bremsstrahlung interpretation requires very large 
Coulomb energy loss rate of ~ 5 x 10 37 ergs s -1 and conse- 
quently the total amount of the energy loss of <~ 10 49 ergs, 
regardless of the gas density and the energy of electrons 
responsible for the hard X-radiation. Given the limited 
energy budget of 7 Cygni, it is more comfortable to at- 
tribute this huge energy release to the acceleration at the 
HXC region, rather than to assume that the electrons are 
supplied by external accelerators. Otherwise, the energy 
requirement would be increased by an order-of-magnitude 
or more, unless we assume very specific arrangement of 
external accelerators around the HXC region. 

If the HXC region contributes only a fraction of the 
EGRET high energy 7-ray flux, the gas density estimate 
of 130 cm~ 3 should be considered as an upper limit. All 
accelerated electrons are presumed to be trapped inside 
the HXC volume. However, it may be possible that high- 
energy (> 100 MeV) electrons effectively escape from the 
HXC, and later interact with very dense clouds resulting 
in high energy 7-rays and synchrotron radio emission, but 
without noticeable X-radiation. The spatial structures of 
the HXC will be comprehensively studied by our forthcom- 
ing Chandra observation. The origin of the 7-ray emis- 
sion from 7 Cygni will eventually be tested by the future 
Gamma-Ray Large Area Space Telescope (GLAST) mis- 
sion. 

We thank Dr. L. Drury for his critical comment concern- 
ing the interpretation of the site of electron acceleration. 
Y.U. thanks Dr. M. Ishida for stimulating discussions and 
his advice on the ASCA data analysis, and also wish to 
thank all the members of the ASCA team. The work of 
Y.U. was supported in part by the Research Fellowships of 
the Japan Society for the Promotion of Science for Young 
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Fig. 1. — (a)-(c): ASCA GIS X-ray images (gray scale) of the 7 Cygni supernova remnant in (a) 0.7—1 keV, (b) 1-3 keV, and (c) 4-7 keV 
energy bands. Two GIS detectors (GIS2, GIS3) are summed. The brightness level indicated at the top of the images is in units of the surface 
brightness of the cosmic X-ray background of the relevant energy interval. Radio contour map (NRAO 4.85 GHz: 7' FWHM) is superposed 
on the ASCA images. A central part of strong emission around an H II region called 7 Cygni nebula at southeast is blanked from the radio 
map. Also shown arc the EGRET 95% confidence error circles (dotted circle: 2EG J2020+4026, thick circle: 3EG J2020+4017). 

Fig. 2. — The GIS radial count rate profile of the clump C2 in the 3—10 keV energy band. A histogram represents the simulated radial 
profile of a point source. Estimated background level is taken from the southern part of the remnant. 

Fig. 3. — A comparison between the energy spectra integrated over the GIS ficld-of-vicw in the 1.2-2.5 keV energy band (top) and the 
3.5-8 keV band (bottom); the south 7 Cygni (filled circle), Field 1 (open circle), Field 2 (rectangle), Field 3 (diamond), and Field 4 (triangle). 

Fig. 4. — The GIS energy spectra extracted from the regions Rl (filled circle) and R3 (open circle), with linear scale. The R3 spectrum is 
multiplied by a factor of 0.3 for display purpose only. The curves show best-fit models, folded through the response function of the instrument. 
The bottom panels plot the residuals of data compared with the thermal emission models. 

Fig. 5. — The GIS energy spectra extracted from the clump C2 (filled circle), where the background data arc taken from the Field 4. The 
curve shows only the thermal component of the best-fit model, folded through the response function of the instrument. The bottom panel plots 
the residuals of data compared with the thermal component. Also shown is the high-energy part of the C2 spectrum whose background data 
are taken from the sourthern portion of 7 Cygni (open circle); the high-energy continuum is hardly affected by the choise of the background 
data. 

Fig. 6. — Broadband spectral power of the 7 Cygni SNR. The range of the power-law fit of the hard X-ray component is shown together 
with 7-ray data (> 100 MeV) of 2EG J2020+4026 taken from Esposito ct al. (1996) and the Whipple TeV upper limit from Buckley et al. 
(1998). The bremsstrahlung photon spectra from the loss-flattened electron distribution are calculated for the electron index s = 2.1 and the 
gas density n = 34 cm~ 3 (long-dashed line), s = 2.3 and n = 130 cm~ 3 (dashed line), and s = 2.3 and n = 10 cm -3 (solid line). The dotted 
lines show the bremsstrahlung spectra corresponding to the acceleration spectra of electrons, i.e. ignoring the Coulomb losses of electrons. 
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Table 1 

Summary of the ASCA archive fields near the 7 Cygni SNR 



Field 


Target 


Coordinate 


Distance a 






(1,6) 


(deg) 


7 Cygni b 


2EG J2020+4026 2 


( 77?92, 2?22 ) 


0.0 


Field 1 


V444 Cyg 


( 76?66, 1?43 ) 


1.5 


Field 2 


NGC 6888 


( 75?55, 2? 42 ) 


2.4 


Field 3 


GRO J2019+37 


( 75? 45, 0?61 ) 


2.9 


Field 4 


GEV 2035+4213 


( 81?22, 1?02 ) 


3.5 



a Angular distance from (l,b) = (77?92, 2?22). 

b This pointing covers the southern part of the 7 Cygni SNR. 
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Table 2 

Results of Spectral Fits to the ASCA data of 7 Cygni 
Parameter R3 R2 Rl C2 CI 



Power Law: 

Photon Index P • • • • • • 1.2 to\ 0.8 ±0.4 1.5 ±0.5 

F 2 -iokev a ••• ••• 1-8 ±8ig 0.98 tlfo 1-7 ±S"i 

Thermal Plasma: 

kT (keV) 0.76 ±^9 0.53 ±0.07 0.56 ±%° 3 s 0.56 flxcd 0.56 fi xod 

EM (10 12 cm- 5 ) b 3.5 to ? 3.0 1.3 t >; ;? 0.47 ±0.02 0.97 ±g"°| 

n c t (10 10 cm- 3 s) 5.8 ±\\l ■■■ 

Neon Line: 

7 Ne (HT 3 P h cm-^- 1 ) c • • • • • • 7.9 + _\% 0.83 ±°;g 5.5 ti;£ 

Photoelectric Absorption: 

AT H (10 22 cm- 2 ) 1.1 ±0.1 1.2 t >; 2 0.84 0.84 flxcd 0.84 flxod 

X 2 1.19(26) 0.80(14) 0.96(34) 1.21(40) 1.25(14) 

a Unabsorbed flux (2-10 keV) in units of 10~ 12 ergs cm~ 2 s~\ 
b Emission Measure: J n c nndV '/AtyD 2 . 

c Photon flux of herium-like Ne K emission line (0.92 keV). 
Note. — Best-fit values and their la errors. 
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(b) 1-3 keV 
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(c) 4-7 keV 



Fig. 1. — GIS images 
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Fig. 2.— Radial Profile of C2 
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Fig. 3. — Background Fields 




Fig. 4-5.— 
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Fig. 6. — Broadband spectrum 



